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ABSTRACT 

The propagation of light of distant sources through a distribution of clumpy 
matter, acting as point mass lenses, produces multiple images that contribute to 
the total brightness of the observed macroimages. In this paper we refine the 
theory of gravitational microlensing for a planar distribution of point masses. In 
the second accompanying paper, we extend the analysis to a three-dimensional 
lens distribution. 

In the two-dimensional case, we derive the probability distribution of macro- 
image magnification, P{A), at high magnification {A — 1 ^ r^) for a low optical 
depth (r ^ 1) lens distribution by modeling the illumination pattern as a super- 
position of the patterns due to individual "point mass plus weak shear" lenses. 
A point mass lens perturbed by weak shear S produces an astroid-shaped caus- 
tic. We show that the magnification cross-section a{A\S) of the point mass plus 
weak shear lens obeys a simple scahng property and provide a useful analytic 
approximation for the cross- sect ion. By convolving this cross-section with the 
probability distribution of the shear due to the neighboring point masses, we 
obtain a caustic-induced feature in P{A) which also exhibits a simple scaling 
property. This feature results in a 20% enhancement in P{A) at A ^ 2/r. 
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In the low magnification {A—1 <S 1) limit, the macroimage consists of a single 
bright primary image and a large number of faint secondary images formed close 
to each of the point masses. The magnifications of the primary and the secondary 
images can be strongly correlated. Taking into account the correlations, we derive 
P{A) for low magnification and find that P{A) has a peak of amplitude ~ 
at A — 1 ^ T^. The low magnification distribution matches smoothly to the high 
magnification distribution in the overlapping regimes A <^ 1/t and A — 1 ^ t^. 

Finally, after a discussion of the correct normalization for P{A), we combine 
the high and low magnification results and obtain a practical semi-analytic ex- 
pression for the macroimage magnification distribution P{A). This semi-analj^ic 
distribution is in qualitative agreement with the results of previous numerical 
simulations, but the latter show stronger caustic-induced features at moderate A 
for r as small as 0.1. We resolve this discrepancy by re-examining the criterion 
for low optical depth. A simple argument shows that the fraction of caustics 
of individual lenses that merge with those of their neighbors is approximately 
1 — exp(— 8r). For r = 0.1, the fraction is surprisingly high: ^ 55%. For the 
purpose of computing P{A) in the manner we did, low optical depth corresponds 
to T < 1/8. 

Subject headings: gravitational lensing 

1. INTRODUCTION 

Gravitational lensing provides a powerful independent tool to probe the distribution of 
matter in the universe and in individual astronomical objects. There are several situations 
in which the lensing objects can be modeled as an ensemble of point masses that produce 
multiple images (microimages) of distant sources. If the angular separations of the microim- 
ages are too small to be resolved observationally, one deals with the macroimage, which is 
the superposition of the microimages. A well-known example of gravitational microlensing 
is the effect of individual stars in a galaxy on the lensing properties of the galaxy (Young 
1981). If the dark matter in the galaxy is in the form of compact objects, the optical depth 
for microlensing could be even higher. Indeed, the microlensing of stars in nearby galaxies is 
being used to probe for the possible existence of massive compact halo objects (MACHOs) 
in our Galaxy (Paczyhski 1986b; Alcock et al. 1993; Aubourg et al. 1993). Another inter- 
esting possibility is that a substantial fraction of the dark matter in the universe is in the 
form of compact objects. Then the universe as a whole has a significant optical depth to 



-3- 



microlensing (Press & Gunn 1973). Finally, microlensing can arise in a globular cluster if 
there are brown dwarfs present (Paczynski 1994). 

The main motivation of this study is to give a systematic theory of gravitational mi- 
crolensing by a random distribution of point masses. There are different regimes of microlens- 
ing and they often require different types of modeling. First, one has to distinguish the cases 
of low (r <^ 1), moderate (r ~ 1), and high (r ^ 1) optical depth, where optical depth r 
is the fraction of the sky covered (locally) by the Einstein circles of all contributing lenses: 
r = T^J^niX^ where rii is the projected (onto the observer plane) surface number density 
of lenses with Einstein radius XE,i (see § 2). Second, the spatial distribution of the lenses is 
usually either "compact" (e.g., stars in a distant galaxy) or "extended" (e.g., cosmologically 
distributed compact objects). A lens distribution is compact if the lenses are distributed in 
a region of scale R which is much smaller than the other distances (those between the source 
and the lenses and between the observer and the lenses) in the problem; it is extended if R 
is comparable to the other distances in the problem. 

In the case of compact lens distributions, the deflection of light rays is generally small 
and essentially occurs in the vicinity of the lens distributions. Consequently, a compact lens 
distribution is usually approximated as a planar (two-dimensional) structure corresponding 
to its projection onto a lens plane that is oriented perpendicular to the line of sight between 
the observer and the source. This "single lens plane" approximation allows us to pose and 
explore well-defined questions concerning the caustic structure of the lens configuration, the 
resulting illumination pattern, the probability distribution of macroimage magnification, the 
temporal variation of the observed brightness of a source (which is the magnification profile 
along a track through the illumination pattern), etc. And while these issues have been 
explored in great detail both analytically and via computer simulations (Kaiser 1992; Mao 
1992; Ranch et al. 1992; see also the book by Schneider, Ehlers, & Falco 1992 and references 
therein), by no means have all the outstanding issues been resolved. 

The greatest progress has been made in understanding the properties of single plane 
lensing under low optical depth conditions. One of the earliest approaches (e.g.. Turner, 
Ostriker, & Gott 1984) advocated modeling the low-r illumination pattern as a superposition 
of what one would obtain for a collection of point masses acting independently. An isolated 
point mass lens produces a circularly symmetric illumination pattern with a divergent spike 
at the origin. Its differential cross-section for magnification A, (Tq{A), is proportional to the 
mass of the lens and scales as A^^' for A ^ 1. A simple superposition of the differential 
cross-sections implies that the probability distribution of macroimage magnification P{A) is 
(Paczynski 1986b) 

P{A) dA = 2TA-^dA for A > 1, (1) 



-4- 



or more generally, 

P{A) dA = 2t{A^ - ly^/'^dA for A - 1 > (2) 

It is well known, however, that neither the r- nor the A-dependence is correct in the 
above distributions, even in the limit of high magnification. As shown by Nityananda & 
Ostriker (1984; also Chang & Refsdal 1984), each lens is subject to an external shear from 
the macroscopic mass distribution as well as its nearest neighbors. An isolated point mass 
lens produces a degenerate point-like caustic. External shear breaks the degeneracy and 
as we shall demonstrate in this paper, the resulting caustic has the shape of an "astroid" 
and its width is proportional to the dimensionless magnitude S of the shear perturbation. 
The corresponding magnification cross-section cr(y4|5') shows a strong feature at A ~ 1/5*. 
This suggests that there ought to be a feature in P{A) at A ~ 1/r since it can readily be 
shown that both the typical magnitude of the macroscopic shear and its random component 
are of order r. The natural extension of the superposition approach, therefore, is to model 
the illumination pattern as a superposition of the patterns due to individual "point mass 
plus weak shear" lenses. Briefly, the resulting probability distribution of magnification is a 
convolution of the magnification cross-section (t(A|S') for a single lens with the probability 
distribution function for the shear p{S). In a subsequent section, we pursue this approach 
and discuss the properties of the P{A) thus acquired. 

Schneider (1987a) adopted a very different approach in attempting to determine the 
asymptotic behavior of P{A) in the limit of high magnification. He argued that the high 
magnification events are dominated by observers lying close to the fold caustics and hence, 
the illumination pattern is dominated by two images of nearly equal brightness. This greatly 
simplifies the problem. Instead of having to identify all the microimages that make up a 
macroimage and then calculate the probability distribution function for their total magni- 
fication, one can simply use the probability distribution function for the individual image 
magnification. The latter can be calculated regardless of the optical depth. This approach, 
however, has one drawback and that is the difficulty in quantifying how large A must be 
before P{A) is well approximated by its asymptotic form. In the low optical depth derivation 
discussed above, the asymptotic form is a good approximation for A ^ t~^. Nonetheless, it 
should be noted that Schneider's derivation is the only rigorous analytic result for finite r. 

Understanding the behavior of P{A) at low magnification (i.e., at 6A = A — 1 ~ r^) is 
also problematic. The distribution quoted above (eq.pf) is not properly normalized as the 
integral / dAP{A) diverges as the lower limit of the integral approaches A = 1. One rather 
cavalier approach is simply to impose a cut-off at some Amin such that the total probability is 
unity. Schneider (1987b) adopted a more sophisticated approach. He argued that a typical 
observer (i.e., one who is not particularly well aligned with any single lens) will see one 
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primary image close to the unperturbed position of the source and a large number of faint 
secondary images (one close to each lens) produced by strongly deflected rays — the so-called 
"diffuse" component. Then the total macroimage magnification associated with the lensing 
event is A = Aq + Yl'h=i where Aq = 1 + SAq is the magnification of the primary image 
and Ai {i > 0) is the magnification of the secondary images. Assuming that 5Ao and Y^iLi 
are statistically independent, Schneider obtained a P{5A) that exhibits a peak at 5A ~ r^, 
is correctly normalized, and has the same A-dependence as equation for A — 1 <^ 1. 
Unfortunately, the two expressions differ in amplitude, a discrepancy that arises because the 
magnifications of the primary and secondary images are not independent. In fact, we will 
show rigorously that the magnifications can be strongly correlated. 

In the above paragraphs, we have hinted at the various calculations and discussion that 
are to be presented in this paper. Our overall goal is to find an expression for P{A) in the 
low optical depth limit that represents a synthesis of the ideas discussed above and that 
is both useful as well as physically motivated. Our approach is a two-pronged one. For 
5 A ^ r^, we perform the superposition of the point mass plus weak shear lenses, allowing 
for the random mutual shear perturbation. In §§ 3.1-2 we calculate the cross-section of an 
individual point mass plus weak shear lens which shows strong caustic-induced features. In 
§ 3.3 we consider the random superposition of the point mass plus weak shear lenses. We 
find that the random shear tends to smear out the strong caustic-induced features of the 
individual lenses and that the deviations from equation around A ~ 1 /r are at the 20% 
level. We also find a simple scaling property for the caustic-induced "bump" for low optical 
depth. For low magnification, we resolve the discrepancy in Schneider's (1987b) calculation 
by allowing for the correlations between the primary and secondary images (§ 4). This gives 
a formula which is valid for A — 1 <^ 1 and which agrees with the expression for 5 A ^ 
in both slope and amplitude. An expression for P{A) that includes all the relevant effects 
is given in § 5.1. We compare this semi-analytic P{A) with the numerical results of Ranch 
et al. (1992) in § 5.2. While the caustic-induced bump in the semi-analytic P{A) agrees 
qualitatively with the bump found in the numerical simulations, the numerical results show 
stronger caustic-induced features at moderate A for r as small as 0.1. The additional features 
are due to caustic configurations more complicated than the simple astroids of the point mass 
plus weak shear lenses. To understand the contribution from the more complicated caustic 
configurations at small r, we consider the leading additional contribution from close pairs 
of point masses with merged caustics. In § 5.3 we calculate the fraction of point masses 
whose caustics are not isolated astroids as a function of r. This gives a simple criterion for 
the optical depth r below which P{A) can be constructed by the superposition of the cross 
sections of individual lenses. 

In the second accompanying paper (Lee et al. 1996; hereafter Paper II), we conduct a 
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similar study for the more complicated situation of a three-dimensional distribution of point 
mass lenses. 



2. LENS EQUATIONS 

We begin by outlining the basic concepts of gravitational lensing, discussing the relevant 
equations and defining our notation. For the latter, we have adopted the notation of Kaiser 
(1992), some of which differs slightly from those in general use. 

For simplicity, we shall only consider gravitational lensing in an Einstein-de Sitter cos- 
mological background weakly perturbed by the gravitational field of the lenses. The line 
element for such a universe is 

cis^ = a^{r^) [(1 + 2(j))dr]'^ - (1 - 2(f)) (dx^ + x^d^f)] ■ (3) 

Let the conformal time ?7 = 1 at the present. The scale factor is a = a^rf, and the Hubble 
parameter is if = Horj'^ with Hq = 2/ao. In this convention, the unit comoving length is 
ao = 6000/i~^ Mpc, where h = ifo/(100kms~^ Mpc~^). The Newtonian (peculiar) poten- 
tial induced by the density inhomogeneities is given by Poisson's equation (in comoving 
coordinates) : 

V20(X) = Qv-'Mx), A(x) = (pix) - P) IP- (4) 

Let us consider an isotropic point source at the origin of our coordinates and an observer 
plane which is perpendicular to the z-axis and at a comoving distance Xso from the source (see 
Fig. 1). In the absence of perturbations, the observer plane would be uniformly illuminated 
and a ray which leaves the source with angle 6 = (^1,^2) would pierce the observer plane 
at cc = Xso^- This defines our planar Lagrangian coordinates x. In the presence of density 
inhomogeneities, the light ray would suffer deflections and the comoving displacement vector 
s{x, x) at distance x from the source is given by 

X 

s{x, X) = -2 / dx' (X - X') V0[a^' + s{x, x'), x'], (5) 


where x' = x'G. For x = Xso, equation (^ defines the mapping from Lagrangian to Eulerian 
coordinates in the observer plane: 

r{x) = x + s{x,Xso)- (6) 



For most compact lens distributions of interest, the deflection of a light ray is small and 
occurs essentially in the vicinity of the lens. In effect, the deflection can be thought of as 
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occurring at a plane that is situated close to the location of the lens(es). In this (single) 
thin-screen approximation, the lens mapping (^) is a gradient mapping with the displacement 
vector s = V\l/, where \1/ is an effective surface gravitational potential. For the particular 
case of single plane lensing by point masses (for simplicity, we assume that all the point 
masses have the same mass m), 

s{x) = -AGm\ ^ = -^i?Z^i 7 

where is the position of the kth. lens projected onto the observer plane (i.e., = Xso^k 
where Ok is the angular position of lens k from the source), is the scale factor at the 
redshift of the lens plane, and xlo, Xsoi and XsL are the comoving distances between the lens 
plane and the observer plane, between the source and the observer plane, and between the 
source and the lens plane, respectively (see Fig. 1). In addition, 

_ / AGrnxLoXsoV'"^ 
xe = (8) 

V O.LXsL J 

is the (Lagrangian) Einstein radius of a point mass on the lens plane. For an ensemble of 
randomly distributed point masses, the optical depth to micro lensing is r = Tmx\, where 
^ = {cllXslIXsoY^/t^ is the projected (onto the observer plane) surface number density of 
point masses and S the physical surface mass density of the lenses. 

The generic behavior and the properties of gradient mappings such as that defined by 
equation (^ are relatively well-known. In fact, the gradient mapping has been used to 
describe a variety of phenomena. The mapping ought to be familiar to the aficionados of 
the Zel'dovich approximation (Zel'dovich 1970) for describing the growth of cosmological 
velocities and density inhomogeneities due to gravitational instability. The evolution of the 
density probability distribution function in the Zel'dovich approximation before significant 
orbit crossing was studied by Kofman et al. (1994). The two-dimensional mapping has been 
used to describe the brightness distribution behind a phase screen (similar to the illumination 
pattern that appears at the bottom of a swimming pool) (Longuet-Higgins 1960; Gurbatov, 
Malakhov, & Saichev 1991). In all cases, however, the results give the brightness distribution 
(or the density of particles) for a single stream. We are not aware of any multiple streaming 
solution for the brightness distribution at the bottom of a swimming pool, for instance. 
The situation in microlensing is quite different in that one must solve the gradient mapping 
equation to locate all the microimages, and then sum the image magnifications to obtain the 
macroimage magnification. This is much more difficult. 



In general, gravitational lensing maps a number of (micro)images with Lagrangian co- 
ordinates a;^, . . . , a;" onto the same Eulerian coordinate r. The observed flux of the image 
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with position a; = a;' is magnified (or amplified) by Ai = l/\V{x^)\, where 

^ _ dr ^ ds 
dx dx 

is the deformation tensor (or Jacobian matrix) associated with the lens mapping and I is 
the identity matrix. If the individual microimages are not resolved observationally, one deals 
with the macro image which has a total magnification (or amplification) A = J2^=i Ai. In 
Lagrangian space {x), the loci of points on which = (i.e., infinite magnification) trace 
out curves that are called critical curves, while the mapping of these curves to Eulerian 
space (r) defines the caustics. In our notation, because we consider a single source and an 
ensemble of observers, the caustics lie on the observer plane. This is different from the usual 
practice of considering a single observer and an ensemble of sources, which has caustics on 
the source plane. 

For a given lens configuration, one can define a differential cross-section a {A) such that 
a{A)dA is the area in the observer plane where the total magnification is between A and 
A + dA. As we shall see, it is also useful to define a "normalized" differential cross-section 
(p{A) = a{A)/ao{A), where cto(A) is the differential cross-section associated with the simple 
lens configuration of a single isolated point mass (see eq.[|lT| below). The function ip{A) is 
a measure of the deviations of the cross-section of a more complex lens configuration from 
that of the single isolated point mass lens. Finally, for lensing by an ensemble of randomly 
distributed lenses, the generalization of the cross-section is the probability distribution of 
macroimage magnification P{A). 



3. CAUSTIC-INDUCED FEATURE AT HIGH MAGNIFICATION 

In § 3.1 we solve the lens equation for a point mass lens perturbed by weak shear. This 
allows us to study analytically the caustic-induced features in the differential cross-section 
(§ 3.2). Finally, in § 3.3, we consider the superposition of the point mass plus weak shear 
lenses and derive the caustic-induced feature in P{A) at high magnification. 



3.1. "Point Mass Plus Weak Shear" Lens: Solution of the Lens Equation 

The lens equation for an isolated point mass lens located at comoving distance xlo from 
the observer and at the origin of the lens plane is 

r = x-xl-^. (10) 
\xr 
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The corresponding differential cross-section is 

= J-%V2 - ^^^^ 

As is well known, the caustic associated with this lens is a degenerate point located at r = 0, 
the degeneracy being due to the spherical symmetry of the lens. 

The degeneracy is lifted if the point mass is not strictly isolated and is perturbed by 
shear, as in the typical lensing configuration in single plane low optical depth microlensing. 
Under low optical depth conditions, lensing (other than low magnification events) is domi- 
nated by the point mass closest to the location where the light ray pierces the lens plane and 
the contribution of the other distant point masses can be treated as perturbations. If we 
choose a coordinate system such that the dominant lens is located at the origin of the lens 
plane and expand the perturbations to the first order in x, the perturbations correspond 
to a constant deflection and a small constant shear and the lens mapping (eqs.0 and 0) 
simplifies to 

r = x + s{x), s{x) ^ dL{x) + a.L + SlX, (12) 

where 

dL{x) = r (13a) 

\x\ 

denotes the influence of the dominant lens, while cxl (the constant deflection) and Sl (the 
shear matrix) represent the perturbative influence of the other lenses on the plane: 



k 



Xk 



Xk 



2' 



(13b) 

In the last equation, S is the magnitude of the shear perturbation and 

f(0,)=f ^"^20. sin 20, \ 

Finally, transforming the Eulerian variable according to r — > r' = r — oti and rotating the 
coordinate axes such that they coincide with the principal shear axes, i.e., 0^ = 0, we can 
cast equation (O) into a more conventional form: 



x-xl-^ + ST{<S)x, (15) 
\xr 
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where 

r(0) = ( ; J ) . (16) 

Equation ([T5|) yields a fourth order equation for the image positions (see Appendix A), which 
has four solutions inside the caustic, and two solutions outside. 

For arbitrary value of S, it is difficult to proceed any further, especially with regards 
to determining the cross-section (t{A) analytically. Fortunately, the magnitude of shear is 
typically S" ~ r and in the limit of low optical depth, is small (see § 3.3). In Figure 2 we show 
the caustic and the contours of constant magnification on the observer plane for a point mass 
perturbed by weak shear {S <^ 1). (The results shown in Figure 2 were obtained numerically 
for 5* = 0.02.) For S* <^ 1, one can observe that the entire four-stream region (i.e., the region 
inside the caustic) has high magnification. The minimum (though high) magnification in 
this region occurs at r' = and it can be found from equations (Al)-(A3) that 

A = (^s - s^y^ ^ s-^ :^ 1 (17) 

at r' = (also Mao 1992). Elsewhere within the caustic, A is even higher and diverges to 
infinity at the caustic. This means that the light rays that delineate the four-stream region 
pass very close to the critical curve in Lagrangian space. To first order in S, the equation 
for the critical line (eq.[A4]) can be expressed as 

a;2 = x|(l + 5cos2?9), (18) 

where we have used the polar representation (x, for the Lagrangian coordinate x. Hence, 
the critical line is a slightly fiattened ellipse. We can study the behavior of the rays that 
pass near the critical line by considering rays with {x/xeY = 1 + i^ and t ^ 1. If we now 
rescale the variables: t' = t/S and g' = g/S = t'/^xeS), equation (Al) reduces at the lowest 
order in 5* to a simple quartic equation for t': 

t'^ - (2 + g'^)t'^ + 2g'H' cos 20 + (1 - g'^) = 0. (19) 

In general, equation ( |19|) has either two or four real solutions (corresponding to the positions 
of the two images outside the caustic or the four images inside). Once the images are 
identified, their magnifications can be summed to get the macroimage magnification factor: 

Therefore, the magnification inside and around the caustic depends on the parameter S in 
the form A oc 1/S. This means that the distribution of A over the observers is a function of 
the combination SA only and any function of A obeys a simple scaling property. 



t'- 



/2 



1) cos'. 



2t' 



t'2 - 2t' COS 20 + 1 



(20) 
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We can obtain a parametric expression for the caustic in Eulerian space by solving 
equations (|19|) and pO|): 



it" - 1) 



t'2 - 2t' cos 20+1 



(21) 



cos 20— sin^ 20 (1 + cos 20) 



1/3 



+ 



sin^ 20(1 - cos 20) 



1/3 



After some tedious algebra, equation ( ^I]) can be simplified to 



f/^ [cos 



2/3 



sm 



2/3, 



22/3_ 



(22) 



Thus the caustic has the shape of an astroid, with the four cusp catastrophes at = (n — 
l)7r/2 [n = 1, . . . ,4) and at a distance g' = 2 (or r' = 2xeS) from the center of the astroid 
connected to each other by four fold catastrophes. The caustic shown in Figure 2, which was 
obtained numerically for 5* = 0.02, is described by equation (521). 



3.2. Analytic Determination of the Caustic-Induced Features in a{A\S) 

One consequence of the scaling relationship noted above is that the caustic-induced 
features in the differential cross-section of a point mass plus weak shear lens scale as 

a{A\S) = aoiAMSA), (23) 

where (Jq{A) is the differential cross-section of an isolated point mass (eq. |jlll) and (p{SA) is 
a function describing the scaling behavior. In Figure 3 we show the "normalized" differential 
cross-section ip{SA) obtained numerically for S = 0.02 (solid line). The modification of 
the cross-section at SA ~ 1 is quite pronounced (a factor of a few higher or lower) and 
the deviations persist over a decade in SA. These features were first found numerically by 
Nityananda & Ostriker (1984; see also Mao 1992). In this subsection, we investigate the 
properties of f{SA) and find analytically the critical values of A and ip{SA). We will also 
provide an analytic fit to ip{SA). 

First it can be seen that ip{SA) tends towards unity asymptotically at high and low mag- 
nifications. As we have noted previously, a{A\S) dA is the area in the Eulerian plane between 
the A and A + dA contours. Essentially, this is equal to the length of the iso-magnification 
contour corresponding to magnification A multiplied by the incremental separation dr' be- 
tween the A and A + dA contours. It is well known that in the vicinity of the fold caustic 
(i.e., in the limit of high magnification A ^ 1/5'), A oc l/i/r^(0) — r', where r^(0) is the 
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distance to the caustic from the center of the astroid. On the other hand, the lengths of the 
iso-magnification contours are nearly independent of magnification because both in size and 
in shape the contours asymptotically approach the astroid-shaped caustic. Consequently, 
o"(74|S') dA oc dr' oc dA/A^. Since the differential cross-section of an isolated point mass also 
depends on A according to (To{A) oc 1/A^ in the high magnification limit, (p{SA) tends to a 
constant value. In fact, to first order in S, the magnification cross-sections of the isolated 
point mass and the point mass plus weak shear cases coincide at high magnification and 

At large distances from the astroid, we can solve equations (Al) and (A2) by only 
considering terms to first order in S and using a perturbation series with respect to 1/r'. 
We find that the resulting iso-magnification contours are given by A ^ 1 + 2{xE/ry + 
2S{xE/r'y cos 2(f). In the limit of low magnification A <^ 1/S (but with A — 1 ^ 5^), the 
contours are similar to those associated with an isolated point mass lens. Consequently, the 
magnification cross-section converges to that of an isolated point mass and ip —>■ 1. 

The interesting features in ip{SA) occur at 5*^4 ~ 1. We can understand these features 
by studying the geometry of the iso-magnification contours on the observer plane (see Fig. 2). 
As we noted above, the iso-magnification contours at large distances from the astroid are 
slightly deformed circles. In going from low to high magnification, the quadrupole moment 
of the contours increases and the iso-magnification contours are continuously deformed. At 
the first critical value Ai, the contour osculates the astroid and due to symmetry, touches 
it at four Eulerian points with coordinates r' = Sxe, <P = {2n — l)7r/4, and n = 1, ... ,4. 
There are four images associated with each of these points, and we can determine their 
Lagrangian coordinates as well as their magnifications from equations (0) and (|20|). Two 
of these images have infinite magnification and the other two have finite magnification; it is 
the latter two that are of interest. From equation (^), we find that the total magnification 
associated with the images of interest is 

Al = 2/3V3S. (24) 

The next contours of A slightly larger than Ai consist of four symmetric arcs around and 
outside the four cusps and because of the loss of the area inside the astroid the corresponding 
(T{A)dA decreases rapidly. This explains the first break in (p{SA) at SAi. 

Next, we consider the iso-magnification contours inside the caustic, i.e., in the four- 
image region. The minimum magnification A2 inside the caustic is located at r' = 0, the 
origin of the Eulerian plane. Using equations ([19|) and (pO]), we find that 



A2 = 1/S. 



(25) 
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[As we mentioned earlier, the exact solution is A2 = 1/5'(1 — S"^) for < < 1, which 
reduces to 1/S for S <^ 1; eq.(]T7|)]. We can also find the magnification in the vicinity of 
the origin by evaluating equations (|TU|) and ( pUD for r' ^ S'x^; or ^' ^ 1. We find that the 



Eulerian position (r', 0) is associated with four images with t'^ 2 = 1 i {r' / Sxe) siiKp and 
= — 1 ± (t'/Sxe) COS0 and that the total magnification is 

Therefore, the differential cross-section is 

r/r' 

aiA\S) = 2Txr' — = 4:7rxlS^e{SA - 1) + 5a{A\S), (27) 

where B is the step function and the term 5a{A\S) ~ A"^/^ (Mao 1992) takes into account 
the contribution to the differential cross-section from contours with A > A2 that are outside 
the caustic and in the vicinity of the cusps (see Fig. 2). The discontinuous jump at A2 



described by equation (|27|) is in good agreement with the jump of height 2 in (p(SA) shown 
in Figure 3. 

Finally, we provide a simple analytic fit (with three pieces for the different A intervals) 
to the numerically obtained differential cross-section of the point mass plus weak shear lens 
shown in Figure 3: 

ip{SA) = 1 + 7.7(5^)3-5 ioiA<Ai, 

- 0-17 , 0.023 A < A < A 

~ {SA - 0.33)1/2 + ^SA - 0.33) tor A^ < A < A2, ^^8) 

= l + %f +(Mf forA>^. 
This fit is shown by the dotted lines in Figure 3. 



3.3. Superposition of Cross-Sections with Distributed Shear 

Let us now consider a collection of point masses (all of the same mass m) randomly 
distributed on a single lens plane in the low optical depth limit (r <^ 1). The combination 
of individual lenses produces several important effects. First, as we discussed previously, 
the lens distribution as a whole induces shear in the neighborhoods of the individual point 
masses. For a point mass perturbed by weak shear, we have found in §§ 3.1-2 that the shape 
of the iso-magnification contours has a universal form — every lens produces a caustic with 
the same astroid shape (Fig. 2) — and that the normalized differential cross-section (f is a 
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function of the combination SA only (Fig. 3). Since the magnitude of the shear S varies 
from lens to lens, there is a distribution of S for the cross-section a{A\S) = (TQ{A)ip{SA) 
near each point mass. In this subsection, we consider the superposition of the cross-sections, 
with the caustic-induced features at different A ~ 1/5*. Other effects of the combination of 
individual lenses will be considered systematically in subsequent sections. 

Under low optical depth condition, the point masses are well separated, and the caustics 
and the high-magnification contours on the observer plane are relatively isolated. Conse- 
quently, the macroimage magnification distribution P{A) at high A can be approximated as 
a superposition of the cross-sections of the individual point mass plus weak shear lenses: 

roc 

P(A) = n dSp(S)a(A\S), (29) 
Jo 

where p{S) is the probability distribution of the shear due to the other lenses on the plane. 
For a random lens distribution with optical depth r = Tinx\^ p{S) is given by (Nityananda 
& Ostriker 1984; Schneider 1987b; Lee & Spergel 1990) 

tS 

P(^^ = (^2^^2)3/2 - (30) 

Note that the shear distribution tends to p{S) = t/S"^ for S ^ t, which is just the shear 
from the nearest neighbor, and that p{S) has a prominent peak at S = t/ \/2 for r ^ 1. 



Substituting equation ( pO]) into equation (p9D, we obtain 

PiA) = nao{A) ^rA) = ^^^^l-^J\(rA), (31) 
where we have introduced the function 

Mr A) = r dy 7?^^ =i+rdy m 



1 + 2/2)3/2 ' « (1 + ^2)3/2 

to describe the caustic-induced feature in the macroimage magnification distribution. Since 
/i(rA) depends on r and A only through the combination tA, the caustic-induced feature 
in P{A) also has a simple scaling property. The function /i(rA), obtained numerically with 
!f{SA) as shown in Figure 3, is plotted in Figure 4 (solid curve). There is a mild "bump" 
located at A ^ 2/t, which is a 20% enhancement at its peak. The bump is similar in shape 
to the features in but it is much weaker and smoother because of the convolution over 
the shear distribution. The function fi{TA) provides a semi-analytic description for the 
caustic-induced bump found numerically by Ranch et al. (1992), but only in the low optical 
depth limit (see § 5). 
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In Figure 4 we also show an analytic fit (dotted curve) to the numerically evaluated 
MrA): 

MrA^-l-OM ^[^^,, (33) 
The maximum deviation of this fit from the numerical result is less than 0.01. 



4. MODIFICATION OF P{A) AT LOW MAGNIFICATION 

The superposition of the cross-sections of individual lenses that results in the caustic- 
induced feature /i(rA) is only valid for the magnifications seen by the observers close to 
one of the astroid-shaped caustics and this, as it turns out, imphes the range A — 1 ^ r^. 
Equation (H) also requires modification at low magnification (i.e., at SA = A — 1 ^ 1) as it 
diverges as 

P(A) ^ (34) 

at y4 = 1 and is not normalizable. 

For the modification of P{A) at low magnification, we are concerned with the observers 
who do not lie particularly near an astroid. These observers see one primary image which 
is barely deflected from the unperturbed position of the source and a large number of faint 
secondary images, one close to each lens, which are produced by strongly deflected rays (the 
diffuse component). If we choose a coordinate system such that a; = for the primary ray, 
then the lens equation near the primary ray is just equation ( |T^ without the cZ^ term, and the 
magnification of the primary image (which we shall denote by Aq) is obtained by summing the 
shear perturbations from the point masses around the primary ray: Aq = {1 — S^)~^ ~ l + S*^, 
where S = \S\ = \J2^k\ and is the shear perturbation from the kth point mass in vector 
notation (see Appendix B). The magnification of the faint image close to the kth point mass 
is Ak = x%/\xk\'^ = Hence the total excess magnification is 

6A = Ao + J2^k~l = \Y.^k\' + Y.\^kf- (35) 

k k k 

Since the lenses have a Poisson distribution, the calculation of P{6A) is well defined. Schnei- 
der (1987b) has obtained an expression for P{6A) under the assumption that the two terms 
in equation (|35D are statistically independent. While that gives a turn-over at 6 A ~ r^, the 
asymptotic behavior at 6 A ^ fails to match the normalization of equation (0). This is 
not altogether surprising since for 6 A ^ r^, both of the sums in equation ( p5D are dominated 
by the nearest lens and the two terms are perfectly correlated. 
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In Appendix B, we determine P{SA) at low magnification, taking into account the 
correlations between the primary and secondary images. This is rather complicated as P{6A) 
formally reads 



P{6A) 



k 



iA-i)-\j:s,\'-Y.\sk\ 

k k 



(36) 



where p{sk) is the probability distribution of the shear from the kih point mass and 5 is the 
Dirac delta function. Although equation ( p6D cannot be reduced to a simple closed form, we 
show in Appendix B that it can be reduced to a single integral, which is a function of the 
combination [A — Vj/r"^ (eq. [B14]). Equation (B14) can be written in the form 



P{A) 



T 



V2{A - 1)3/2 



91 



{A -I] 



2r 



(A2 - 1)3/2 



9i 



{A 



(37) 



where we have used the fact that A+1 for 5 A <^ 1. The correction function gi [{A — 1)/ t'^] 
takes into account the effect of the combination of individual lenses for small magnification. 
Numerical integration of equation (B14) gives the function gi and the distribution P{A) 
plotted in Figure 5. An excellent analytic fit to gi is given by 



9i{y) 



ye 



-Tr/iy 



(7r/23/2 + y)- 



(3^ 



This fit has the same asymptotes (eq. [B15]) as the integral (B14), and it does not differ from 
the numerical integration results by more than 0.6%. 

The distribution (^) has the following properties. It is a normalized distribution: 
/ dAP{A) = 1. It has a sharp peak P{A) ^ 0.16/r2 at 6 A ^ 0.84r2 (see Fig. 5). For 
A — 1 ^ T^, it has the asymptote 



PiA) 



V^(A- 1)3/2' 



(39) 



which matches the low magnification {6 A <^ 1) asymptote (eq.[|^) of the distribution (0). 
Since both distributions (PT| ) and ( P7| ) give the same result in two overlapping regimes: 
A -C 1/r and A — 1 ^ r2, it is easy to construct a macro image magnification distribution 
that takes into account the effects of both the caustics and the diffuse component: 



PiA) 



2t 



f,{rA)g,[{A-l)/r' 



(A2 - 1)3/2 

In the appropriate limits, this distribution becomes either equation (PT| ) or ( P7| ) 



(40) 
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5. FINAL FORM OF P{A) AND COMPARISON WITH 
NUMERICAL SIMULATIONS 

5.1. Renormalization and the Final Form of P{A) 

In § 3.3 we assumed that the surface density of the astroid-shaped caustics on the 
observer plane is the same as the projected surface density of point masses (see eq. |2^). There 



is, however, a net convergence of the light rays by the overall lens distribution. This increases 
the density of astroids by a factor A, where A = {1 — r)~^ is the average magnification on 
the observer plane. Therefore, there should be an additional factor A in the normalization 



of P{A) in equation (^). In the very low optical depth limit (r <^ 1), the correction due to 
this factor A ^ 1 + 2r is small. However, for optical depths of order 0.1, the correction is 
not negligible. 

As we noted in § 1, Schneider (1987a) derived an analytic result for P{A) in the high 
magnification limit, which is exact even for finite optical depths r. The starting point of 
his analysis is to consider the probability distribution of the shear on the lens plane p{S) 
(eq.|]30|). Since the magnification for a light ray that encounters shear S on the lens plane 



is y4 = |1 — 5^1^^ and each light ray corresponds to an image in Lagrangian space, the 
distribution p{S) can be transformed into the probability distribution function, Pl{A), for 
the magnification of the individual images in Lagrangian space. For high magnification, 
S = 1 + 6S with \6S\ = 1/2A, so 



Pl(A)=p{S = 1] 



d6S 



dA 



2(1 + ^2)3/2 



A-\ (41) 



To obtain the probability distribution function for the magnification in Eulerian space P{A), 
we have to first multiply Pl{A) by the factor A/ A, which takes into account the focusing 
of the light rays by the lenses (Schneider 1987b; Lee & Spergel 1990), and then double the 
magnification [A 2A), which takes into account that the high magnification events are 
dominated by a pair of very bright microimages that form close to a fold catastrophe. The 
final result is 

^^^^ " (1 + r2)3/2 ' = (1 _ r)2(l + r2)3/2 ^ (^2) 

which has the usual A~^ form. Although this is an exact analytic result for finite r, it 
is limited in applicability to sufficiently high magnification such that the macroimage is 
dominated by two microimages. As we saw in § 3, for low r, one actually has to go to 
fairly high magnification {A ^ 1/t) for this to be a good approximation. The normalization 
in equation (^) differs from the normalization 2r in equation (|40| ) by the factors A and 



^1 + r2)-3/2. For r ^ 0.1, A ^ 1 + 2t and (1 + t'^)~^/^ ^ 1 - 3x^2, and the correction to 
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the normalization 2r is mainly due to the factor A. Note that we have already derived the 
factor A with the argument given in the previous paragraph. 

We can now collect together the three effects of the combination of individual lenses 
discussed above — the caustic-induced feature from equation (0), the low magnification 
modification due to the diffuse component (eq. |3^), and the renormalization due to the net 
convergence of light rays and multiple streaming (eq.^^) — and find that for r ^ 1, the final 
form of the macroimage magnification distribution due to a two-dimensional distribution of 
point masses is 

PiA) = (,_^),',"^^,)3/2 - ^r'^' hi-A) g, [{A - l)/r^] . (43) 



5.2. Comparison with Numerical Simulations 

Several authors (e.g., Paczyhski 1986a; Ranch et al. 1992; Wambsganss 1992) have 
used numerical simulations to calculate the macroimage magnification distribution, P{A), 
produced by a two-dimensional distribution of point masses. In particular, distributions with 
high resolution in A have been obtained from Monte Carlo simulations by Ranch et al. for 
r = 0.1, 0.2, and 0.3. We now compare the semi-analytic P{A) derived in this paper, equation 



(^), to these numerical results. As noted, our derivation assumes low optical depth (i.e., 
r ^ 1). It is not immediately obvious whether optical depths in the range 0.1 < r < 0.3 are 
sufficiently small and whether we should expect a good agreement between the semi-analytic 
and numerical results. In fact, as we shall see, r = 0.1 is not sufficiently small. There 
are differences between the semi-analytic and numerical results that are indicative of finite 
optical depth effects. Ideally, comparison should also be made for smaller r, but numerical 
P{A) with the required accuracy is not available [the computational requirement increases 
rapidly with decreasing r because the caustic-induced feature shifts to higher A (oc 1/r) 
while the amplitude of P{A) at large A (oc r) decreases]. 

In Figure 6 we show the high resolution P{A) obtained by Ranch et al. for r = 0.1 and 
0.2 (histograms). The data are those shown in Figure 3 of Ranch et al. (1992), but they are 
plotted in the form {A^ — lY^'^P{A) and the bin sizes at large A are slightly different. For 
comparison, the solid lines show the full semi-analytic P{A) (eq. ^3|), and the dotted lines 
show the semi-analytic P{A) with the caustic-induced feature only (i.e., eq. 03] without the 



function gi). Note that in a plot of {A^ — l)^/^P(y4), the semi-analytic distributions without 
the low magnification modification gi or the caustic-induced feature /i are simply horizontal 
lines of amplitude 2r(l — 'r)~^(l + r^)~^/^ (dashed lines). 



At logy4 < 0.1 (or 6 A < 0.3), the semi-analytic and numerical results for r = 0.1 are 
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in good agreement, but the results for r = 0.2 are slightly different in shape and amplitude. 

Since this range of A is dominated by the low magnification modification gi. wc conclude that 
the function gi derived in § 4 is valid for r < 0.2. (Recall that the peak in P{A) produced 
by gi is located at 6A ^ 0.84r^, which is much less than 0.3 for both cases.) The semi- 
analytic caustic-induced feature also provides a reasonably good fit to the numerical results 
at log A > 1 and 1.2 for r = 0.1 and 0.2, respectively. At intermediate A, there are significant 
differences between the numerical and semi-analytic results. The numerical distribution is 
lower than the semi-analytic distribution in the neighborhood of the minimum at A 2.5, 
and it is higher than the semi-analytic distribution at smaller and larger A. This pattern of 
deviations is similar to the cross-section shown in Figure 3 and suggests that the function 
fi derived in § 3 underestimates the strength of the caustic-induced feature at intermediate 
A. It is, however, important to note that the differences between the semi-analytic and 
numerical results decrease in strength with decreasing r and should be relatively small for 
r ^ 0.05. In the next subsection we shall analyze the cause of the differences at intermediate 
A. 



5.3. Re-examining the Criterion for Low Optical Depth 

In our analysis of the caustic-induced feature in P{A) for low optical depth (§ 3), we 
assume that all the point masses are well separated. Then the dcffcction near each point mass 
is due to that point mass and the shear perturbation from the other lenses, and an astroid- 
shaped caustic is associated with each point mass. However, since the lenses are randomly 
distributed on the lens plane, some of the lenses have close neighbor (s), separated by less 
than a few Einstein radius xe- In these cases, the caustics are not isolated astroids but more 
complicated structures produced by the collective effect of two (or more) point masses. If 
the surface density (or optical depth) of the lens distribution is low enough that the fraction 
of lenses with close neighbor(s) is small, the contribution to P{A) by these configurations is 
negligible, and the analysis in § 3 is valid. For larger r, however, we have to take into account 
the more complex configurations and can do so by evaluating the macroimage magnification 
distribution as a series: P{A) = Pi{A) + P2{A) -\ — •, where Pi{A) is the contribution by the 
point masses perturbed by shear, P2{A) is the contribution by close pairs of point masses, 
etc. The first term Pi{A) is the distribution derived in § 3 (but with a maximum cutoff in 
the convolution over shear). The second term P2{A) can be evaluated (approximately) as a 
convolution of the cross-section a{A\d) for two point masses separated by distance d with 
the probability distribution for d if we consider a point mass and its nearest neighbor as a 
two-point-mass lens and neglect the perturbation from the other lenses. 
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To understand the contribution to P{A) from close pairs of point masses, we must first 
look at some of the properties of lensing by two equal point masses on a single lens plane (see 
§ 3 of Paper II and Schneider & Weiss 1986 for details). For consistency with the notation in 
Paper II, we shall express the Lagrangian separation d between the lenses in units of ^/2xe- 
In the limit c? S> 2, the region near each of the point masses is perturbed by the weak 
shear {S = l/2dP) from the other point mass, and there are two astroid-shaped caustics. 
The caustics move towards each other (and become asymmetric) as d decreases, and they 
touch when d = 2. As d decreases below 2, the number of caustics changes from two to one 
to finally three (see Fig. 3 of Paper II). In Paper II, this sequence of caustic topologies is 
denoted as topology types A' (for d>2), B' (for 1/^2 < d<2), and C (for d < 1/^2). In 
all cases, the normalized differential cross-sections ^p{A) are qualitatively similar to that for 
the point mass plus weak shear lens, but the caustic-induced features can be significantly 
stronger (compare Fig. 4 of Paper II to Fig. 3 of this paper). As in the point mass plus weak 
shear case, there is a discontinuous jump in (p{A) at the minimum magnification Amin inside 
the caustic(s). In the limit d 3> 2, the minimum magnification inside the two astroid-shaped 
caustics is A^i^{d) ^ A2{S = l/2d^) = 2d'^ (eq.[|5l). However, as Witt & Mao (1995) have 
shown, Amin is a non-monotonic function of d and has a global minimum of 3 when d = \f2. 
If we now consider the convolution of these cross-sections with the probability distribution 
for c?, it is clear that the resulting distribution P^i^Ai) should show an enhancement in the 
caustic-induced feature near A = 3. 

In their analysis of the caustic-induced feature in -P(v4), Ranch et al. (1992) compared 
the results from the full scale Monte Carlo simulations to the results from a simpler two- 
point-mass model (see their Fig. 6 for the r = 0.2 case). Their two-point-mass calculation 
is in fact an approximate numerical evaluation of -Pi (A) + P2(^)- They found that the two- 
point-mass model is able to reproduce partly the caustic-induced feature at intermediate A. 
In particular, the distribution P(A) from the two-point-mass model shows a dip at A 3. 
Ranch et al. suggested that the slightly lower value {A ^ 2.5) of the dip location found in 
the full Monte Carlo simulations is due to non-negligible contributions from configurations 
of three or more point masses. 

Is there a simple explanation for the relatively strong contribution to P{A) by the 
collective effect of two (or more) point masses for r as small as 0.1? For a Poisson distribution 
of point masses, the probability that the nearest neighbor to a point mass is at a distance 
less than d (again in units of ^/2xe) is simply 1 — exp(— 27rnx|;(i^). If we ignore the deflection 
due to the other lenses, the point mass and its nearest neighbor is a two-point-mass lens. 
As we mentioned earlier, a two-point-mass lens produces two astroid-shaped caustics only if 
d > dAB, where d^B = 2 is the separation at which the caustic topology changes from type 
A' to B'. Therefore, a simple estimate for the fraction of point masses whose caustics are 
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not isolated astroids is 

Pna(r) = 1 - exp(-27™x|4^) = 1 - exp(-8r). (44) 

For the contribution to P{A) from P2{A) (and higher order terms) to be neghgible, Pna{T) ^ 
1 or r ^ 1/8. It is immediately clear that r = 0.1 is not sufficiently small for P2{A) to be 
negligible. For r = 0.1, we estimate that about half of the point masses produce caustic 
structures that are more complex than the astroid shape: Pna = 0.55. This is consistent with 
the illumination pattern shown in Figure 1 of Ranch et al. (1992). In § 5.2, we concluded from 
the comparison that the contribution from P2iA) should be reasonably small for r <^ 0.05; 
this corresponds to Pna ^ 0.33. 

6. SUMMARY 

In this paper we have attempted to build upon the various approaches developed in 
previous studies in order to develop a systematic theory of gravitational microlensing by a 
planar distribution of point masses in the low optical depth limit. In particular, we have 
derived a practical semi-analytic expression for the probability distribution of macroimage 
magnification, P{A). 

At high magnification {6 A ^ r^), we model the illumination pattern as a superposition 
of the patterns due to individual "point mass plus weak shear" lenses. The shear perturbation 
near each point mass is induced by the neighboring point masses. It breaks the degeneracy 
of the caustic of an isolated point mass lens and produces an astroid-shaped caustic. The 
convolution of the magnification cross-section of the point mass plus weak shear lens with 
the probability distribution of shear yields P{A) = 2t{A^ — l)^'^/^/i(ryl), where the function 
fiijA) (eq.[^; Fig. 4) describes the caustic- induced feature in the macroimage magnification 
distribution. In effect, fi{TA) introduces a mild "bump," a 20% enhancement, aX Ak. 2/t. 
Since fiirA) depends on r and A only through the combination rA, the caustic-induced 
feature in P{A) exhibits a simple scaling property. To facilitate future computations, we 
have provided a useful analytic fit to fiirA) (eg. p3[|). 

We should point out that the results derived in §§ 3.1-2 for the point mass plus weak 
shear lens may also have applications in the analysis of gravitational microlensing by physical 
binary systems. Microlensing searches towards the Galactic bulge and the Large Magellanic 
Cloud have already discovered microlensing events by close binaries with a single merged 
caustic (Udalski et al. 1994; Alard, Mao, & Guibert 1995). Since wide binaries are more 
common than close ones, there should be a significant number of events due to wide binaries 
(Di Stefano & Mao 1996). If the components of the binary are sufficiently far apart, the 
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region near each point mass is perturbed by the weak shear from the other point mass. In 
these cases, there are two separate astroid-shaped caustics, and the results derived in §§ 3.1-2 
are apphcable. 

At low magnification {5 A <^ 1), the macroimage consists of a bright primary image 
barely defiected from the unperturbed position of the source and a large number of faint 
secondary images (the diffuse component) formed close to each of the lenses. The magni- 
fications of the primary and the secondary images can be strongly correlated. Taking into 
account the correlations, we find that P{A) = 2r(A^ — — l)/r^], where gi (eqs.pTj 
and [B14]; Fig. 5) represents the low magnification correction. An excellent analytic fit to 
Qi is given by equation (^). The function gi prevents P{A) from diverging at A = 1, and it 
introduces a sharp peak of amplitude P{A) ^ 0.16/r^ at A — 1 ^ 0.84r^. 

We also discussed the renormalization of P{A) due to the net convergence of light rays 
and multiple streaming. Finally, collecting together the above results, we find that in the 
low optical depth limit (r <C 1), 



In order to determine the realm of validity of the above semi-analytic expression, we com- 
pared it against P{A) obtained from Monte Carlo simulations by Ranch et al. (1992) for 
r = 0.1 and 0.2. At low magnifications (log A < 0.1), we find that the r = 0.1 semi-analytic 
and numerical results are in good agreement with each other. At greater optical depths, 
differences arise both in shape and amplitude. We, therefore, conclude that the low mag- 
nification modification gi is valid for r < 0.2. The numerical and the semi- analytic results 
are also in good agreement in the high magnification regime (log A > 1 and 1.2 for r = 0.1 
and 0.2, respectively). At intermediate A, however, the semi-analytic expression does not 
match the numerical result even for optical depth r = 0.1. The deviations arise because 
the function /i derived in § 3 underestimates the strength of the caustic-induced feature 
at intermediate A. The deviations tend to diminish with decreasing optical depth, and we 
expect them to be relatively small for r <^ 0.05. 

In order to understand the discrepancy between the numerical and semi-analytic results 
for r as small as 0.1, we re-examined our derivation of the function fi describing the caustic- 
induced feature in P{A). In our derivation, we assume that a unique astroid-shaped caustic 
is associated with each point mass. In a random distribution of point masses, there may 
arise groups of two (or more) lenses that lie sufficiently close to each other and give rise to 
caustic configurations more complicated than isolated astroids. If the surface density of point 
masses is small, such lens configurations will be rare. If, however, the surface density is large, 
they will be more common and their contributions cannot be neglected. For a distribution of 
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optical depth r, we estimate that the fraction of point masses whose caustics are not simple 
astroids is Pna.{T) — ^ — exp(— 8t). For r — 0.1, about half of the point masses fall into this 
category. For the purpose of computing the macroimage magnification distribution P{A) in 
the manner we did, a simple criterion for low optical depth is Pna -Clorr-Cl/S. 

In the accompanying Paper II, we extend our analysis to the more complicated situation 
of a three-dimensional distribution of point mass lenses. 
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A. LENS EQUATIONS FOR POINT MASS PLUS ARBITRARY SHEAR 



Expressing r' and x in polar coordinates (r', 0) and (x, d) respectively, equation (|T5|) 
can be written as a fourth order equation for the Lagrangian coordinate x in terms of the 
Eulerian coordinate r'\ 

{C-lf[C-{2 + Q^)C + l]= -2Sg' cos 20 C'(C - 1) + S'g\' + 2S\\C -if- S\\ (Al) 



tawd 



tan( 



(A2) 



where we have defined ( = {x/xe)^ and g = r'/xE- In this notation, the Jacobian of the 
mapping is 



dr' 



dx 



1 25 cos 2^9 . 

'-e — c — ^■ 



(A3) 



The critical line (i.e., the loci of Lagrangian points where the Jacobian \V\ vanishes) is given 
by 

S 



C 



cos 2^9+ Wcos2 2?9 + 



52 



(A4) 



which corresponds to the so-called "Cassini oval." The associated caustic is non-trivial and 
it is not a point singularity. Based on equations (Al) and (A2), we expect four images in 
the interior of the caustic and two outside. 



B. LOW MAGNIFICATION LIMIT OF P{A) 

Let us consider a random distribution of = imR^ point mass lenses within a circle of 
Lagrangian radius i?, where n is the projected surface number density of lenses. The total 
shear at the origin due to the N point masses is given in matrix form by equation ( |13b| ). 
We shall find it convenient to write the two distinct components of the shear matrix Sl as a 
vector: S = J2k=i ^k, where Sk = {xE/\xk\Y (cos20fc, sin20fc) is the shear perturbation from 
the kth point mass. 

To determine the low magnification limit of P{A), we need to evaluate the probability 
distribution function for the following combination of shears (eq. |P5||): 

N N 

5A = A-l = \Y.Skf + Y.\skf. (Bl) 

fc=l k=l 

The first term in equation (Bl) is just S"^ = jSp, and the probability distribution of this 
term alone is found easily from the distribution (pO]). The second term in equation (Bl) is the 
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so-called diffuse component, and its probability distribution has previously been considered 
by Schneider (1987b). The evaluation of the probability distribution of the sum of these two 
terms, which is 6A, is more complicated. 

Let us start with the formal expression for the probability distribution of 6 A as a function 
of the statistically independent shears (eq.[p6|): 



P{SA) 



N 



N 



k=l 



k=l 



(B2) 



where p{sk) is the probability distribution of the shear from the kth point mass and 5 is 
the Dirac delta function. Equation (B2) can be simplified if we express each term in the 
integrand as a Fourier transform. For p{sk), we have 



d^tkq{tk) e 



itk-Sk 



where the characteristic function q{tk) is (Nityananda & Ostriker 1984) 



r 



q{tk) = 1 - —4, 



(B3) 



(B4) 



and r = Txnx\ is the optical depth. Because of the cross terms in \S\^, it is inconvenient 
to use the direct Fourier transform of the delta function in equation (B2). Instead, one can 
introduce an auxiliary integration to express the delta function in the form 



5 



N N 

{A-i)-\Y.Sk?-Y.\^k\ 

k=l k=l 



I 



d?X6 
1 



1X1 



N 1 / ^ \ 

k=l J \ k=l J 



(B5) 



(2vr)= 



J d?X J dv J (fu exp 



N 



N 



-IV[W 



lU- 



[x-y: sk] 



k=l 



k=l 



where 5^ is the two-dimensional delta function and w = A — 1. Substituting equations (B3) 
and (B5) into equation (B2), we can perform the integration with respect to Sk- 



/ (fsk exp[ivsl - i(tk - u)-si,] = — exp{-i\tk - u\'^/4:v), 
J V 



and then the integration with respect to tk'. 
1 



I dHkq{tk) (^^) exp(-2|tfc - u\yAv) = 1 



N 



(B6) 



(B7) 
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where 



f{v,u) 



(B8) 



ITT 



and Jo is the zero-order Bessel function. In equation (B8), z = iu^/Sv, 

g{z)^[{l + 2z)Io{z) + 2zh{z)]e-\ 

and Iq and Ii are the modified Bessel functions. 

We can now group together the N integrations over tk (eq.[B7]) to obtain 



(B9) 



N 



-f{v,u) 



(BIO) 



in the hmit N ^ oo. Then 

1 



P{SA) 



J (fx JdvJ (fu exp [-iv{w - X^) - iu-X - f{v, u)] . (Bll) 



In equation (Bll), the integration with respect to X is 

J V 



(B12) 



and we are left with two integrations with respect to v and u. Substituting equations (B8) 
and (B12) into equation (Bll), and using variables v and z instead of v and u, we obtain 



2 roo roo 

P(5A) = - dze-^' dvcos 

TT Jo Jo 



WV — T 



exp 



(B13) 



Finally, after performing the integration with respect to we have P{5A) in the form 



P{5A) 



w^/^ Jo 



dzg{z) exp 



T^T 9/ N 



(B14) 



where w — 5A and g{z) is defined in equation (B9). Equation (B14) must be integrated 
numerically in general, but it has the following analytic asymptotes: 



P{5A) 



2^ y^l/2^-n/iy 



y/2T 



1 -U-3/2 
2 y 



for y = 6A/t'^ < 1, 
for y = 6A/t^ > 1. 



(B15) 
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Fig. 1. — A schematic diagram illustrating our notation for single plane gravitational lensing. 
A light ray that leaves the source with angle 6 has Lagrangian coordinates x. It is deflected 
on the lens plane and reaches the observer plane at Eulerian coordinates r. 

Fig. 2. — The caustic and the iso-magnification contours on the observer plane for the point 
mass plus weak shear lens. They were obtained numerically for shear S = 0.02. The axes 
are in units of xeS. The caustic has the shape of an astroid with the four cusps at a distance 
2xeS from the origin. The contours of critical values Ai = 2/3^/3S and ^42 = 1/5' are also 
indicated. 

Fig. 3. — The "normalized" differential cross-section, ip{SA) = (T{A\S)/ao{A), of the point 
mass plus weak shear lens {solid line). It was obtained numerically for shear S = 0.02. The 
dotted lines show the analytic fit, equation (pSf) . 

Fig. 4. — The function fi^rA) (eq. |^) that describes the caustic- induced feature in the 
macroimage magnification distribution, P{A), at low optical depth r. The solid line is the 
numerical result using (p{SA) shown in Fig. 3, and the dotted line is the analytic fit, equation 

Fig. 5. — (a) The function gi[y = {A — 1)/t^] (eq.[0) that describes the modification of 
P{A) at low magnification in the low optical depth limit. (6) The macroimage magnification 
distribution P{A) (= gi[y]/ V2y^/^T^) at low magnification in the low optical depth limit. 

Fig. 6. — Comparison of the semi-analytic macroimage magnification distribution, P{A), 
with the Monte Carlo results of Ranch et al. (1992) for optical depth r = 0.1 and 0.2 
{histograms). The solid lines are the semi-analytic P{A) with the caustic-induced feature 
and the modification at low A (eq. |^), and the dotted lines are the semi-analytic P{A) with 



the caustic-induced feature only (i.e., eq. |^ without the function gi). 
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